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Abstract. High-resolution IRAS maps are used to search for 
the presence of stellar-wind bow-shocks around high-mass X- 
ray binaries (HMXBs). Their high space velocities, recently 
confirmed with Hipparcos observations, combined with their 
strong stellar winds should result in the formation of wind 
bow-shocks. Except for the already known bow-shock around 


Vela X-1 (Kaper et al. 1997), we do not find convincing evi¬ 


dence for a bow-shock around any of the other HMXBs. Also 
in the case of (supposedly single) OB-runaway stars, only a mi- 
norit y appe ars to be associated with a bow-shock (Van Buren 
et al. 1995). 


We investigate why wind bow-shocks are not detected for the 
majority of these OB-runaway systems: is this due to the IRAS 
sensitivity, the system’s space velocity, the stellar-wind prop¬ 
erties, or the height above the galactic plane? It turns out that 
none of these suggested causes can explain the low detection 
rate (~40 %). We propose that the conditions of the inter¬ 
stellar medium mainly determine whether a wind bow-shock 
is formed or not. In hot, tenuous media (like inside galactic 
superbubbles) the sound speed is high (~100 km s“^), such 
that many runaways move at subsonic velocity through a low- 
density medium, thus preventing the formation of an observ¬ 
able bow-shock. Superbubbles are expected (and observed) 
around OB associations, where the OB-runaway stars were 
once born. Turning the argument around, we use the absence 
(or presence) of wind bow-shocks around OB runaways to 
probe the physical conditions of the interstellar medium in the 
solar neighbourhood. 
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1. Introduction 

OB-runaway stars are massive, OB-type stars with high space 
velocities. This small (~20 % of the O stars, less than 5 % of 
the early B stars), but significant excess of high-velocity ob¬ 
jects suggests a syste matical prod uction of fast movers (Feast 
& Shuttleworth 1965 , Stone 1979 ). Among the OB runaways, 
velocities well above 100 km s ~^ ha ve been o bserve d (Gies 
Conlon et al. 


bility to remain bound after the supernova explosion of the ini¬ 


& Bolton 1986 


199C 


sive star (e.g. Van Oijen 1987). 


The two most popular scenarios for the production of run¬ 
aways seem to operate at roughly the same rate (Hooger- 
werf et al. 2001) and are both capable of producing run¬ 
aways with velocities up to 200 km s“^ (Portegies-Zwart 200C 


Leonard 1995); (i) The binary supernova scenario (Blaauw 


1961, Zwicky 


tially most massive star (Boersma 1961). If it remains bound. 


Blaauw |1993[ ), which 
is roughly ten times the average space velocity of “normal” 
OB stars in the Milky Way. Due to their high velocities, OB- 
runaways may have travelled great distances since their for¬ 
mation. The criteria used to classify OB stars as “runaways” 
are (Blaauw 1961): (i) A high space velocity; often a threshold 
space velocity of 30 km s“^ is used, roughly corresponding 
to three times the typical space velocity of OB-type stars, (ii) 
The reconstructed path of the runaway should start in a “par¬ 
ent” OB association. Since the definition by Blaauw, the term 
OB runaway has often been used to include all of the high- 
velocity OB stars, irrespective of the identification of a parent 
OB-association. A large distance above the galactic plane can 
also be used as an indication for the runaway-nature of a mas- 


the binary, now consisting of a massive, rejuvenated main- 
sequence star and a compact remnant (a neutron star or a black 
hole), will travel through space with a high velocity. As soon as 
the massive star becomes a supergiant, material from its dense 
stellar wind (or through Roche-lobe overflow) is intercepted by 
the gravitational field of the compact star and accretes, power¬ 
ing a strong X-ray source; a high-mass X-ray binary (HMXB). 
This scenario thus predicts that all HMXBs are runaway ob¬ 
jects. Recently, this prediction has been co nfirm ed with Hip- 
parcos observations (Chevalier & Ilovaisky 1997, Kaper et al. 


1999). This also suggests that several supposedly single OB 


runaways might have a, so far undetected, compact compan¬ 
ion. For a review on HMXBs and OB-runaway stars, see Kaper 

(|200lD. 


The high space velocity and powerful stellar wind of the 
OB-runaway star will give rise to a strong interaction with the 
ambient medium. In the case of supersonic movement, a wind 


bow-shock is formed (Baranov et al. 1971). The wind bow 


1957), where a massive star in a binary receives 


a high space velocity after the supernova explosion of its (ini¬ 
tially) more massive comp anion; (ii) The cluster ejection sce¬ 
nario (Poveda et al. 1967[ ), where dynamical interactions in a 
compact cluster result in the ejection of one or more of its mem¬ 
bers. Especially encounters between binary systems are ef fec- 
tive in producing high-velocity objects (e.g. Mikkola 1983 ). 

The current version of the binary supernova scenario in¬ 
ches a phase of mass transfer (Van den Heuvel & Heise 


1972): when the initially most massive star fills its Roche lobe 


(e.g. when becoming a supergiant) mass is transfered to its 
companion which eventually will become the most massive star 
in the system. As a consequence, the system has a high proba- 
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shock accumulates gas and dust from the interstellar medium 
which is heated by the ultraviolet radiation field of the OB star. 
The heated dus t subs equently radiates at infrared wavelengths 
(Draine & Lee 1985 ) and the shocked gas becomes visible in 
strong optical e mission lines like Ha and [O III] 4959,5007 A 
(e.g. Kaper et al. 1997 ). Data obtained with the InfraRed Astro¬ 
nomical Satellite (IRAS) has been used to search for the pres¬ 
ence of wind b ow-sh ocks around OB runa ways b y Van Bu- 
ren & McCray ( 1988 ) and Van Buren et al. ( 1995 ). They de¬ 
tected excess infrared (60 pm) emission for about 30 % of the 
188 OB stars in their sample. For a significant fraction promi¬ 
nent arc-like features could be resolved using high-resolution 
maps (Noriega-Crespo et al. 1997). When spatially resolved, 
the shape of the bow-shock indicates the direction of motion of 
the system. The stand-off distance of the shock provides a con¬ 
straint on the space velocity, the stellar-wind parameters, and 
the density of the ambient medium. 

In this study we investigate whether wind bow-shocks are 
as common for HMXBs as for (single) OB-runaway stars, by 
searching for their infrared emission in high-resolution maps 
produced from data obtained with the IRAS satellite. We find 
that also in the case of HMXBs only a minority form a wind 
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Table 1. Prediction of the angular separation between wind bow-shock and HMXB. The HMXB peculiar velocity (i.e. corrected 
for differential galactic rotation and peculiar solar motion) in the radial and tangential direction is listed in column 6 and 7, 
respectively; the resulting space velocity Ugpace is given in column 8 . The first five columns list the identification of the X-ray 
source, the name and spectral type of the OB companion (if known), the distance r, and the height z above the galactic plane. The 
predicted (and observed, final column) angular distance between OB star and wind bow-shock are given in the last columns. For 
the HMXBs that have not been observed with Hipparcos (ftan)pec = 50 km s“^ has been adopted. If a question mark is added, 
infrared emission is detected, but no clear arclike structure could be resolved. The prediction of the stand-off distance is based on 
stellar-wind parameters taken from Puls et al. (1996) and Kaper et al. (1998); for Be stars it was assumed that M = 10“®Mq yr“^ 
and Vac = 250 km s“^. 


Object 


Spectral Type 

r 

(kpc) 

Z 

(pc) 

('^rad)pec 

(kms“^) 

('t^tan )pec 

(km s"^) 

'?^space 

(kms-i) 

jpred 

^shock 

(pc) 

rpred 

^shock 

(') 

bow-shock? dshock 

(') 

0114-H650 


B0.5 Ib 

3.8 

171 

-9.6 

29.4 

31.0 

2.6 

2.4 

no 


1223-624 

GX301-2 

B1.5 la 

5.0 

-3.0 

-1.6 

50 

50 

2.8 

1.9 

7 

1.5(7) 

1700-377 

HD 153919 

06.5 laf 

1.7 

65 

-42.0 

57.3 

71.0 

3.1 

6.3 

7 

8.0(7) 

1907-H097 


B0.5 Ib 

4.0 

37 


50 

50 

0.6 

0.5 

7 

0.5(7) 

1538-522 

QV Nor 

BO lab 

5.5 

206 

-81.9 

50 

96 

1.9 

1.2 

no 


1119-603 

Cen X-3 

06.5 II-III 

8.0 

47 

16.3 

50 

53 

1.7 

0.7 

no 


1956-H350 

Cyg X-1 

09.5 lab 

2.5 

135 

-10.9 

41.4 

42.8 

5.0 

6.9 

no 


0900-403 

Vela X-1 

B0.5 Ib 

1.8 

123 

-22.3 

38.5 

44.5 

1.4 

2.7 

yes 

1.0 

0236-H610 

V615 Cas 

BOe 

2.0 

42 

-30.2 

22.5 

2,1.1 

.14 

.24 

no 


0535-H262 

V725 Tau 

09.7e II 

2.0 

-90 

-41.5 

57.5 

70.9 

.07 

.12 

no 


0352-H309 

XPer 

09e III-IV 

0.8 

-235 

-32.0 

13.4 

34.7 

.51 

2.2 

no 



bow-shock. We investigate whether the detection of a bow- 
shock relates to (i) the kinematical properties of the runaways; 
(ii) the height above the galactic plane; or (iii) the spectral 
type of the OB runaway. The remaining option, namely that the 
physical conditions of the ambient medium mainly determine 
whether a wind bow-shock is formed, seems to be the most 
likely one. If true, this provides the opportunity to use OB run¬ 
aways as probes of the physical conditions of the interstellar 
medium. 


2. A search for bow-shocks around HMXBs 


In our search for the presence of wind bow-shocks around 
high-mass X-ray binaries we selected a number of systems for 
which the conditions for the formation of a bow-shock seem 
favourable: a large space velocity and an OB companion with a 
strong stellar wind and a high luminosity. In practice this means 
that our sample consists mainly of systems hosting an OB su¬ 
pergiant (cf. Kaper |1998|). We also included a few Be/X-ray 


binaries (with early-type OBe companion). The HMXB sam¬ 
ple is listed in Table To determine the size of the area on the 
sky to inspect for the presence of a bow-shock, we estimated 
the angular distance from the star at which a bow-shock should 
form. The stand-off distance (dg) is set by the balance in ram 
pressur e betw een the stellar wind and the ambient medium (cf. 
Wilkin 1996 ). Using typical values for the wind mass-loss rate 
M (in Mq yr“^), wind terminal velocity v^o and space veloc¬ 


ity u* (in km s ^), and the density (in cm ^) of the ambient 
medium n(= p/p,mH), the stand-off distance can be written 


ds — 


0.9pc 


(M/10 ®M 0 /yr) x (uoo/10^km/s) 


( 1 ) 


(t;*/50km/s) V n/cm~^ 

We calculated dg for the systems in our sample (Table |^. 
Stellar-wind parameters were taken from Puls et al. ( 1996[ ) and 
Kaper et al. (1998). The space velocity is obtained from the ob¬ 
served Hipparcos proper motion and radial velocity, corrected 
for differential galactic rotation and peculiar solar motion (cf. 
Moffat et al. 1998| ). The interstellar medium density was com¬ 
puted with the mean density model of Dickey & Lockman 
( 199C ). The typical area on the sky that we inspected for the 
presence of a bow-shock is 0.5° x 0.5°. 


2.7. HIRAS maps of the HMXB environment 


We used the program HIRAS (Bontekoe et al. 1994) to produce 


high-resolution IRAS maps of the circumstellar environment of 
the HMXB sample, to inspect the area for the presence of wind 
bow-shocks. The IRAS satellite scanned almost the whole in¬ 
frared sky in four wavelength bands (12, 25, 60, and 100 pm) 
in strips of about 6' wide. The HIRAS software takes advan¬ 
tage of the fact that a given part of the sky has been observed 
several times, resulting in strips that partly overlap. Combin¬ 
ing these observations then yields a significant enhancement in 
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Fig. 1. HIRAS image of the circumstellar environment of 
Vela X-1 at 60 jjm. The arc-like structure of a wind bow-shock 
is clearly resolved. The arrow indicates the direction of the sys¬ 
tem’s space velocity. East is to the left and North is up; the im¬ 
age size is 0.5° x 0.5°. Infrared emission of the bow-shock is 
also detected in the other three IRAS bands. 


Fig. 2. HIRAS image of the environment of Cygnus X-1 at 
60 /rm. The feature located roughly 20' east of Cyg X-1 is part 
of a larger filament that outlines the Cygnus X region (Cash 
et al. 1980); a large X-ray emitting region possibly heated by 
numerous supernovae and winds from massive stars. East is to 
the left and North is up, the image size is 1.0° x 1.0°. 


spatial resolution. Eor details of the HIRAS-algorithm and the 
used method we refer to Bontekoe et al. (1994). The resolution 
of the images depends on the local sampling pattern and, there¬ 
fore, varies with position on the sky; in our HIRAS images it 
typically is in the range from 0.4' to 1.7'. 

Eor each system HIRAS images were prepared in the four 
IRAS bands. A wind bow-shock has been detected in a narrow- 
band Ha image of the field surrounding the B-supergiant 
HD77581, the companion of the X-ray pulsar Vela X-1 (Kaper 
et al. 1997). The bow-shock is clearly present in the HIRAS im¬ 
ages too. Eigure shows a HIRAS image of Vela X-1 at 60 pm, 
with the direction of its peculiar motion indicated. We prepared 
images in the other three IRAS bands (12, 25 and 100 pm, not 
shown here), but the bow-shock (and at 60 and 100 pm the H II 
region surrounding the system) is detected in all four bands. 
The arc-like structure of the bow-shock is resolved at 25 and 
60 pm. Images obtained with ISOCAM onboard the Infrared 
Space Observatory confirm the detection at 6 and 12 pm, and 
reveal the fine structure present in the Ha image (Kaper et al. 


1998 ). The symmetry axis of the bow-shock is parallel to the 
system’s direction of motion; its space velocity is 45 km s“^ 


(based on Hipparcos observations; note that Kaper et al. 1997 
wrongly quote 90 km s“^). Moreover, the infrared emission 
of the shockfront is coincident with the region producing Ha 
emission. 

Based on the size of the H II region surrounding HD77581, 
Kaper et al. (1997) derive that the number density of the ambi¬ 
ent medium is ~10 cm“^, which is about ten times high er than 
predicted with the model of Dickey & Lockman ( 1990 ) and a 
distance of 120 pc above the galactic plane. The higher density 



Fig. 3. HIRAS image of the area on the sky centered at 1700- 
377 at 60 pm. Many infrared sources are located in the vicin¬ 
ity of 1700-377, making it difficult to identify a possible wind 
bow-shock feature associated with this HMXB. East is to the 
left and North is up, the image size is 0.5° x 0.5°. 


is also consistent with the observed stand-off distance, which 
is less than predicted (Table [^. 

None of the fields we inspected show evidence for arc-like 
structures that can be associated with a wind bow-shock around 
any of the other HMXBs. It is not surprising that the three 
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Fig. 4. HIRAS image of GX301-2 at 60 ^m. Infrared emission 
is detected which could be attributable to a wind bow-shock. 
However, this is a complicated area on the sky, many other in¬ 
frared sources are located within the field of view. East is to the 
left and North is up, the image has a size of 0.5° x 0.5°. 



Fig. 5. HIRAS image of the HMXB 1907H-097 at 60 pm. Possi¬ 
bly this system is producing a wind bow-shock, but this could 
be a chance superposition. Infrared emission was detected in 
all four IRAS bands. Note that the proper motion of 1907H-097 
is not known. East is to the left and North is up, the image has 
a size of 0.5° x 0.5°. 


Be/X-ray binaries, which do not drive strong stellar winds, do 
not produce a bow-shock. On average, these sys tems a lso have 
lower space velocities (Chevalier & Ilovaisky 1997 ); this is 
consistent with the predictions of models describing the evo¬ 


lution of massive binary systems (Van den Heuvel et al. [2000 ). 
However, a few OB-supergiant systems do show some interest¬ 
ing features that are worthwhile discussing. Eigures ^ - 1| show 
HIRAS images of some of these systems. The size (and direc¬ 
tion) of the depicted arrow is proportional to the peculiar space 
velocity of the system (the size of the arrow corresponds to the 
displacement over the next 60,000 years). 

In three cases it remains unclear whether the observed 
HMXBs are producing bow-shocks. In the case of 1700-377 
and GX301-2 bright objects are present in the field possibly 
hampering the identification of a bow-shock (Eigs. || and ^ 
60 pm images). GX301-2 is smTounded by a faint ring-like 
structure at an angular distance of 2' towards the north. As 
we do not know the proper motion of the system it is hard to 
say whether this structure is related to GX301-2. 1907H-097 
shows some emission in all four IRAS-bands at the location 
of the source, a bow-shock might be present. At a distance of 
^ 4 kpc the shockfront is expected to have an angular separa¬ 
tion of ^ 0.5', consistent with what is observed. 

In summary, only one of the HMXBs in our sample is, be¬ 
yond any doubt, associated with a wind bow-shock. The ques¬ 
tion remains why the other systems do not produce a (de¬ 
tectable) bow-shock? They are moving at high velocity and 
have strong stellar winds powering an X-ray source. Unfortu¬ 
nately, our sample is too small to perform a statistical study 
to determine why wind bow-shocks are a bsent. The sample of 
OB runaways studied by Van Buren et al. ( 1995 ) is much better 
suited for such an investigation. 


3. Wind bow-shocks around supposedly single 
OB-runaway stars 


Van Buren et al. (1995) conducted an extended search for wind 


bow-shocks around OB runaways. Based on pie-Hipparcos 
proper-motion and radial-velocity measurements a total of 188 
candidate runaway stars was compiled from the catalogs of 
Garmany et al. ( 1982|), Cruz-Gonzales et al. (1974) and the en¬ 
tire list by Stone ( 1979| ). Van Buren et al. (1995) correctly point 
out that this list should not be considered systematic or com¬ 
plete since it inherits the biases of the input catalogs. Among 
that sample, 58 candidate bow-shocks were identified based on 
low-resolution images from the IRAS Sky Survey Atlas (ISSA). 


In a follow-up study by Noriega-Crespo et al. (1997) the 58 


candidate bow-shocks were more closely examined. They an¬ 
alyzed 60 pm IRAS high-resolution (HiRes) images of these 
objects that were produced with a maximum correlation algo¬ 


rithm (Eowler & Aumann 1994), a method similar to the one 


we discussed in the previous section. Objects were classified as 
being either (i) bow-shock producing objects, or, (ii) if a shock- 
front detection was inconlusive but the object showed excess 
infrared radiation, as excess IR objects. Three objects did not 
fall in any of these catagories and were rejected (or rather, clas¬ 
sified as non-detections). 


As mentioned above. Van Buren et al. (1995) concluded 


that about a third of the 188 OB runaways are associated with 
a bow-shock, i.e. the majority of the systems do not produce an 
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(observable) wind bow-shock. Our aim is to determine whether 
the properties of OB runaways with a bow-shock differ from 
those for which a bow-shock seems to be absent. We investigate 
the spatial distribution and the kinematical properties of the OB 
runaways, their galactic height distribution, and their spectral 
types. 


3.1. Properties of Van Buren’s OB-runaway sample 


We checked the Hipparcos database and found that 168 stars 
from Van Buren’s total of 188 candidate runaways are included. 
We omitted stars with spectral type other than O or B and re¬ 
jected all objects located at distances beyond 4 kpc because of 
the large uncertainty in proper motion. As a result we were left 
with a sample of 148 candidate runaways, to which will be ref- 
ered to as “Van Buren’s sample”. Radial velocities were taken 
from the SIMBAD database; for objects for which no radial ve¬ 
locity is measured, the tangential velocity is multiplied by a 
factor y^3/2 when calculating the space velocity. Proper mo¬ 
tions and radial velocities were corrected for differential galac¬ 
tic rotation and the peculiar solar motion. We followed the ap¬ 
proach described in Moffat et al. ( 1998| ) which uses a flat galac¬ 
tic rotation curve. 


3.1.1. Spatial distribution 

The spatial distribution of the OB runaways included in 
Van Buren’s sample (and with a distance less than about 4 kpc) 
is shown in Fig. The top panel displays the OB runaways 
with a detected wind bow-shock (filled triangles) as well as 
those exhibiting excess infrared emission (filled circles). The 
density contours were calculated using the binormal kernel 
method, where the position of each runaway star is replaced 
by a binormal distribution (i.e. with standard deviation in 
both X and y). Subsequently, the convolution of all binormal 
distributions results in a continuous OB-runaway density distri¬ 
bution in the galactic plane. For instance, let di be the distance 
between an arbitrary position {xi, yt) in the galactic plane and 
the location {x, y) of an OB runaway: 

di = \/{x - Xi)'^ + [y - yiY ■, (2) 

then the kernel is: 

1 --^4 

K,{d,/a,) = —e ^ (3) 

ZTT 

where is an adaptive spread in the kernel function. We chose 
the value of (the width of the distribution) to depend on the 
heliocentric distance (r) to the respective runaway: 

CTr = a X r (4) 

We varied a between 0.1 and 0.5, resulting in a more contin¬ 
uous distribution with increasing a; we adopted a = 0.3 to 
produce Fig. Even though the positional accuracy of each 
runaway is considerably smaller along the tangential direction, 
(Tr is taken to be symmetric in both directions for reasons of 
simplicity. Adding all individual kernels, the density distribu¬ 
tion becomes: 

n ^ 

Pix, 2/) = X! —Kr{d^/(Jr) (5) 

1 ^ Y 

1—1 




O no IR excess 
• IR excess 
A bowshock 


Fig. 6. The spatial distribution of the OB runaways in Van Bu¬ 
ren’s sample (bottom panel). The Sun is located at the origin 
of the diagram. The density contours were calculated with the 
so-called kernel method: each location of a runaway is replaced 
by a normal distribution with a standard deviation correspond¬ 
ing to 30% of the heliocentric distance. The spatial distribution 
of the OB-runaway stars that produce either bow-shocks (filled 
triangles) or show excess infrared emission (filled circles) is 
displayed in the top panel. The space density of runaways with 
a bow-shock decreases significantly with distance, most likely 
due to the reduced probability of detecting a bow-shock. 


Note that each term carries a normalization factor with the di¬ 
mensions kpc~^. The distribution of p{x,y) thus reflects the 
surface density of OB runaways. 

The same method was used to derive the density distribu¬ 
tion of bow-shock producing runaways only (top panel Fig. 

It was assumed that objects that show excess infrared emis¬ 
sion are producing wind bow-shocks. Fig. ^ indicates that the 
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number of OB runaways that produce observable bow-shocks 
drops significantly beyond a distance of 2 kpc. The relative 
abundance of objects showing either a bow-shock or IR-excess 
remains nearly constant 40 %) up to a distance of about 
2 kpc. Within 0.5, 1.0, and 2.0 kpc we respectively find 13 (out 
of 33), 18 (44), and 39 (98) runaways with a bow-shock; this 
cotTesponds to 39,41, and 40 %, respectively. Between 2.0 and 
4.0 kpc we encounter a bow-shock for 9 out of 50 runaways 
(18 %). Therefore, even if the sample becomes less complete 
at larger distances, the roughly constant fraction of bow-shocks 
among the runaways up to a distance of 2 kpc suggests that the 
sampling of the runaway population producing bow-shocks be¬ 
comes incomplete beyond that distance. This might be due to 
the reduction in the detection sensitivity of IRAS. 

A short intermezzo: The sensitivity of IRAS is, at best, 
~ 0.2 Jy at 12, 25, and 60 ^m, and 1 Jy at 100 /rm, 
per O'.25 x O'.25 or O'.5 x O'.5 pixel at 12 and 25 or 60 
and 100 /rm, respectively (Moshir et al. 1989). The highest 
IRAS flux per pixel measured in the bow-shock associated 
with Vela X-1 (Fig. |) is 0.03 (0.01), 0.41 (0.01), 2.59 (0.06), 
and 1.04 (0.21) Jy/px at 12, 25, 60, and 100 /rm, respectively; 
the value between brackets is the level of the infrared back¬ 
ground in that direction. The heliocentric distance to Vela X-1 
is 1.8 kpc; according to the quoted IRAS sensitivity the bow- 
shock would still be marginally detectable (only at 60 /im) 
if the distance to Vela X-1 is increased by about a factor 3. 
In terms of integrated flux (2.8, 38, 95, and 34 Jy) the bow- 
shock of Vela X-1 is relatively bright. For comparison, the 
integrated flux of the “text-book-example” bow-shock around 
a Cam (distance 1.1 kpc) is 0.8, 45, 310, and 470 Jy (Van Bu- 
ren et al. 19951 The angular extent of the bow-shock around 
Vela X-1 is ~ 8'; at a 3 times larger distance it would just be 
resolvable in a HIRES image. Given the strong variations in 
the galactic infrared background emission, it is expected that 
several bow-shocks have remained undetected in the IRAS sur¬ 
vey. However, bow-shocks as bright as those around Vela X-1 
and a Cam should be detectable up to several kpc. 

A striking feature of Fig. ^ is the anisotropy of bow-shock 
producing objects. A high concentration of bow-shocks and IR- 
excess objects is found within 2 kpc from the Sun in quadrant II 
(90° < I < 180°), while in the opposite quadrant IV (270° < 
I < 360°) the lowest fraction of bow-shocks is detected. The 
dumbell shape apparent in the contours of the bottom panel of 
Fig. ^ reflects the spiral-arm structure of the Milky Way. 


3.1.2. Space velocities 

In Fig. ^ the space velocity and distance z with respect to the 
galactic plane are plotted of the OB runaways in our sample. 
This figure shows that most bow-shocks are detected around 
objects with velocities below 40 km s“^. Of the 76 objects that 
have a space velocity less than 40 km s“^, 14 are associated 

* The latter values indicate that the dust in the bow-shock around 
Vela X-1 is hotter than that around a Cam. The X-ray source could 
provide additional heating. 


76 



(Vsp.ce) [km/s] 


Fig. 7. Space velocity versus distance z to the galactic plane. 
This graph shows that only a small number of objects with a 
very high space velocity or with a large z distance produce 
bow-shocks (the meaning of the symbols is indicated in the 
legenda). On the other hand, a large fraction of the objects 
with a space velocity less than 40 km s“^ do show a bow- 
shock: 42 % are associated with a bow-shock or show IR ex¬ 
cess; at higher velocities this rate is only 22%. The number of 
objects within a given space-velocity or z-distance bin is indi¬ 
cated along the axes. The cumulative number of objects with a 
bow-shock (black), with IR excess (grey), or none of the two 
(white) is included. 


with a bow-shock and 18 show an IR-excess, together making 
up for about 42% of the total, while among the 72 objects that 
are faster than 40 km s“^ a lower fraction is found (6 bow- 
shocks and 10 objects with IR-excess, resulting in a 22% de¬ 
tection rate). At space velocities above 80 km s“^ none of the 
objects could be associated with a wind bow-shock, and only 
three (out of 18) were observed with excess infrared emission. 
Apparently, a very high space velocity is not a favourable con¬ 
dition for the production of a wind bow-shock. This observa¬ 
tion is in agreement with the prediction on the basis of numer¬ 
ical simulations of bow-shocks by Comeron & Kaper ( 1998| ) 
that at very high velocity a shockfront becomes unstable and 
disrupts. 

Runaways with large separations from the galactic plane 
also show a low bow-shock detection rate. Among the objects 
that have a z distance larger than 250 pc only 2 out of 13 were 
observed with an IR-excess (15% detection rate). This low de¬ 
tection rate of bow-shocks at large galactic heights can be un- 
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Fig. 8. Hertzspmng-Russell diagram of Van Buren’s sample of 
runaways. Only objects located within 250 pc from the galactic 
plane and with a space velocity in the range of 5 to 80 km s“^ 
are included. There is no obvious correlation between the de¬ 
tection of a bow-shock and the spectral type of the OB-runaway 
star. 


derstood in terms of dust content. Since the dust content of the 
interstellar medium decreases with galactic height (Van Steen- 
berg & Shull 1988), a runaway too far from the galactic plane 
may not show an observable bow-shock because not enough 
dust is being swept up. 

Both the space velocity and the height above the galactic 
plane are physical parameters that have to be taken into ac¬ 
count when considering whether a bow-shock is observable or 
not. However, the majority of runaways without bow-shocks 
are not located at extreme distances from the galactic plane, 
nor have very high space velocities (Fig. Obviously, other 
factors must play an important role as well. 


3.1.3. Dependence on spectral type 


Besides the space velocity of the OB-runaway star and the 
density of the ambient medium, the stellar-wind parameters 
M and Voo determine the physical conditions under which a 
bow-shock is formed (cf. Eq. 1). The wind mass-loss rate M 
is strongly dependent on the luminosity (Mcx Earners & 


Leitherer 1993, Puls et al. 199C). i e. supergiants have much 


stronger winds than main sequence stars. Furthermore, the visi¬ 
bility of the wind bow-shock (at infrared wavelengths) depends 
on the irradiation of the swept-up dust by the ultraviolet flux 


from the massive star. One would expect that supergiants are 
more capable of forming a detectable bow-shock in compari¬ 
son to main sequence stars. 

In Fig. ^ the Hertzsprung-Russell diagram of the OB- 
runaway sample is presented; the symbols indicate whether a 
bow-shock (or IR excess) has been detected. The OB-runaway 
sample contains a relatively large fraction of supergiants and 
giants. This might be due to a selection effect; Hipparcos ob¬ 
served only the brightest stars on the sky. Inspection of the dia¬ 
gram learns that the detection rate of bow-shocks is about 35 %, 
independent of the luminosity class of the OB runaways: 

- From a total of 36 supergiants (luminosity class I, la, Ib, 
lab and II), six objects show bow-shocks and seven show 
excess infrared emission (36 % detection rate). 

- The sample contains a total of 44 giants (luminosity class 
III and IV). Among these, seven show bow-shocks, and 
seven objects show excess infrared emission (32 % detec¬ 
tion rate). 

- Out of 38 main sequence stars (luminosity class V), there 
are six with bow-shocks and eight with infrared excess 
(37 % detection rate). 

Also the effective temperature of the star does not seem to 
affect the detectability of a wind bow-shock. 


4. On the absence of wind bow-shocks around runaway 
objects 

Although we found some observational indications that the 
presence of a wind bow-shock depends on the space velocity of 
the OB runaway and its height above the galactic plane, these 
hndings cannot explain why the majority of systems do not pro¬ 
duce a wind bow-shock. In the following we explore the pos¬ 
sibility that perhaps some runaways are moving subsonically, 
despite their high space velocity. With space velocities found 
of several tens of km s“^ this would require extreme condi¬ 
tions of the interstellar medium, conditions that are known to 
be present in hot bubbles and galactic coronal gas. 

The (isothermal) speed of sound in the interstellar medium 
(ISM) has a value of ^1 km s“^ at a temperature of 100 K, 
^10 km s“^ at 10"* K and ~100 km s“^ at 10® K. Runaway 
stars thus move supersonically unless they travel through re¬ 
gions with temperatures of the order of 10® K. Such hot re¬ 
gions are present in the Milky Way, the so-called hot bubbles, 
and appear e.g. as holes in the distribution of neutral inter¬ 
stellar matter (Heiles |l979|jl984|[I^ Shull & Saken |1995|). 


The hot regions are likely produced by the combined action of 
stellar winds and supernovae in a central star cluster (i.e. an 
OB association), which sweeps the surrounding gas and dust 
away and confines this material into dense filaments (walls). 
The ambient density of the medium is thereby reduced to val¬ 
ues ~ 10“® cm“®. The temperature may reach 10® K, result¬ 
ing in the production of diffuse X-ray emission, such as ob¬ 


served from the Cygnus superbubble (Cash et al. 1980). The 


ionized gas inside the hot bubble can also be observed at radio 
wavelengths, with a synchrotron component originating near 
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0 no IR excess 
• IR excess 
^ bowshock 


0 OB association 
^ HI shells, superbubbles 


O no IR excess 
• IR excess 
^ bowshock 


0 OB ossoclotlon 

0 HI shells, bubbles (Heiles 1979, 1984) 
^ New superbubble (Heiles 1998) 


Fig. 9. The distribution of H I shells in the galactic plane, in the 
range 10<1<250 and -10<b<10 according to Heiles ( 1979 , 


1987) (shaded regions). The fourth quadrant of the galactic 


plane has not been included in these studies. The Sun is lo¬ 
cated in the origin of the diagram; the direction of the galac¬ 
tic cente r is in dicated. OB associations (taken from Melnik & 
Efremov 1995| ) are represented by circles including a plus sign 
(see legenda). OB runaways and their paths (the path length 
corresponds to the projected distance travelled during the last 
10 Myr). The large concentration of OB associations in the 
fourth quadrant suggests that also here hot bubbles must be 
present. Most of the OB runaways producing a bow-shock do 
not seem to be contained inside a hot bubble. 


the walls. The Sun is located in the local hot bubble, a hot 
medium with a deficiency of neutral matter. The walls of this 
bubble are f airly w ell defined by the distribution of interstellar 
dust (Lucke 1978 ) an d the detection of interst ellar a bsorption 
lines (Frisch & York 1983 ). Recently, Heiles ( 1998 ) reported 
the discovery of a nearby superbubble which may have merged 
with the local bubble, forming a highly oval cavity elongated 
towards 1~230. 

To test the hypothesis that such hot interstellar gas indeed 
prevents runaways from producing wind bow-shocks, we in¬ 
vestigated the distribution of runaways in relation to the distri¬ 
bution of hot bubbles in the galactic plane. We used the H I sur¬ 


veys carried out by Heiles (1979, 1984), and compiled a list of 


H I shells, supershells and other shell-like objects. The surveys 
span a range in galactic longitude 10° < I < 250° and latitude 
-10° < b < 10°. The recently discovered nearby superbub¬ 


ble (Heiles 1998), next to the Local Bubble, is also included. 


Fig. 10. Zoom-in on the central area of Fig. y. Some well- 
known hot bubbles are labelled. 


The studies by Heiles do not cover the fourth quadrant of the 
galaxy (270° < I < 360°). The large number of OB associa¬ 
tions in that area make it very likely that hot bubbles are present 
there as well. At longitudes 317° and 327° filaments extending 
over several degrees are dete cted i n the Australian Mononglo 
(MOST) survey (Green et al. 1999 ). 

Figures ^ and show the locations and trajectories of run¬ 
aways (the length of the track corresponds to the distance trav¬ 
elled in 10 Myr), projected onto the galactic plane, with re¬ 
spect to H I shells and bubbles. Figure E (a zoom in on the 
central region in Figure focusses on the solar neighbour¬ 
hood (within 1.3 kpc). OB associations are represented as cir¬ 
cles with a cross in the middle (see lege nda); i n Fig. 1^ the 
associations comprised in De Zeeuw et al. ( 1999 ) are included. 
Apart from Heiles’ New Superbubble and the Eridanus Loop, 
the shape of the H I shells is assumed to be spherically symmet¬ 
ric. The estimated distance and angular size in galactic latitude 
and longitude are used to determine the radius of the bubble. 
A large degree of uncertainty exists in the size of the bubbles 
along the line of sight. 

Although both the OB-runaway and the hot-bubble distri¬ 
bution are not accurately known, the OB runaways produc¬ 
ing bow-shocks are mostly found in regions devoid of hot 
bubbles. Remarkable is the asymmetry in the distribution of 
OB runaways producing bow-shocks with respect to the Sun: 
quadrant II (in the direction of Cepheus) contains significantly 
more bow-shocks than the other quadrants. This anisotropy 
agrees with the distribution of hot bubbles. In the nearby (su- 
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0 no IR excess 0 OB association 

• IR excess ^ HI shells, superbubbles 

^ bowshock 

Fig. 11. The distribution of the HMXBs in our sample is shown 
with respect to hot bubbles and OB associations. Each HMXB 
shows a track that reflects the distance travelled in 10 Myr; in 
some cases these tracks identify the parent OB association of 
the system. Cyg X-1 and X Per are very likely enclosed in hot 
superbubbles. Also 1700-377 and V615 Cas may be situated 
in such hot regions, although the presence of bubbles at the 
positions of 1700-377 and V615 Cas is uncertain. 


per)bubbles we hardly encounter any OB runaway with a bow- 
shock. 

Dove & Shull (1994) estimate that the average volume Ail¬ 
ing factor of (super)bubbles and cavities in the galactic plane is 
50-80%, which is consistent with the fraction of bow-shock 
producing objects among the runaway stars in Van Buren’s 
sample. Apparently, the soundspeed in these regions is high 
enough to prevent the formation of wind bow-shocks. 


4.1. The HMXB sample 


Can the detection of a wind bow-shock around only one HMXB 
(out of 11) also be explained by the distribution of hot bub¬ 
bles? Figure |ll| displays the distribution of the HMXBs in our 
sample. Cyg X-1 is located close to the edge of the Cygnus 
superbubble, while Vela X-1, the only HMXB for which a 
bow-shock clearly has been detected, seems to move through a 
medium devoid of hot bubbles. At least four HMXBs are likely 
to be travelling through a hot bubble. 


1700-377: Benaglia & Cappa (1999) report the presence of 


an H 1 structure that may be due to a wind-blown bubble created 
by the O-supergiant companion of 1700-377, HD 153919. The 


H I region increases in size towards higher galactic latitudes, 
forming some kind of horseshoe shape. The mean velocity of 
the bubble corresponds to a distance ~2 kpc, in agreement with 
the distance of 1700-377 (cf. Ankay et al. 2001). 

V615 Cas. According to Figure this HMXB is located 
near the edge of a large superbubble; two other bubbles are 
found at roughly the same longitude. Shell GS 135H-08-34 has 
a larger catalogued distance than V615 Cas (3.1 and 2.0 kpc, 
respectively), but taking into account the uncertainty in these 
distance determinations, they could possibly be associated with 
each other. On the other hand, the center of GS 135H-08-34 lies 
at a larger angular separation from the galactic plane than V615 
Cas (+8° as compared to +1°). It is not quite clear whether 
the extent of GS 135H-08-34 is sufficient to enclose V615 Cas. 
Heiles (19841 estimates an angular size of only ^ 6° of GS 
135H-08-34 in galactic latitude. 

X Per. Two H l shells can be discerned in the vicinity of 
X Per. The smaller of the two bubbles (GS 165-21-4) is cen¬ 
tered around galactic latitude b= —21° corresponding to a sep¬ 
aration from the galactic plane of ~ 260 pc at a distance of 
0.72 kpc (Heiles 1984). The extent of this shell along the 2 - 
direction is estimated to be ~160 pc at that distance. X Per is 
situated 235 pc below the galactic plane, and is therefore likely 
enclosed in this bubble. 


Cyg X-1. This HMXB is located in the boundary region of 
the Cygnus superbubble (see Figure Cash et al. (1980) per¬ 
formed observations with HEAO-1 of this region and showed 
that the observed X-rays emitted from this region is outlined 
by a ring of elongated Ha filaments, suggesting a common 
origin. Moreover, interstellar absorption considerations yield a 
distance of about 2 kpc of the X-ray feature, thus at roughly the 
same distance as Cyg X-1. The temperature of the X-ray emit¬ 
ting plasma contained in the Cygnus superbubble is estimated 
to be 2x10® K. 


5. Discussion 

One of the main conclusions of our work is that the formation 
of a wind bow-shock around an OB-runaway star strongly de¬ 
pends on the temperature and density of the ambient medium. 
Apparently, the majority of the OB runaways are moving 
through hot bubbles distributed in the galactic plane. These hot 
bubbles are created by the combined action of stellar winds and 
supernovae of OB stars, most of them members of OB associ¬ 
ations. Assuming that all OB runaways originate in OB asso¬ 
ciations, it is very likely that they travel through such hot bub¬ 
bles for a substantial fraction of their (relatively short) lifetime. 
In the following we use the absence of a wind bow-shock as 
a diagnostic to probe the distribution of hot bubbles in the so¬ 
lar neighbourhood. Finally, we present arguments why HMXBs 
should have a higher probability to be located inside a hot bub¬ 
ble than dynamically ejected OB runaways. 
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Fig. 12. The fraction of bow-shock producing runaways (“bow- 
shock probability density”) is shown together with the observed 
distribution of local hot bubbles. The bow-shock probability 
distribution is determined with objects that are located within 
300 pc from the galactic plane and have space velocities be¬ 
tween 5 and 80 km s“^. The HMXB sample is also included, 
yielding a total sample of 137 runaways. In the top panel, the 
dotted contour follows a runaway number density of 1 kpc“^, 
marking the region outside of which the bow-shock probability 
density is determined by less than 1 object per square kpc. 


5.7. The distribution of hot bubbles in the solar 
neighbourhood 

Various people have tried to map out the extent of the local bub¬ 


ble (see Frisch 1995 and references therein). Due to its high 


temperature, some nearby runaways that are moving through 
this bubble are not producing wind bow-shocks, even though 
they have substantial space velocities. The presence, or rather 
the absence, of bow-shocks can thus serve as a diagnostic to 
constrain the extent of the local bubble. An attempt to map 
out hot galactic regions based on the relative presence of bow- 
shock producing runaways is shown in Fig. E- The bow-shock 
probability density reflects the fraction of runaways producing 
a bow-shock. The galactic distribution is derived using the bi- 


Fig.l3. As the bottom panel of Fig. [^, but now calculated 
adopting = 0.15r to illustrate that the obtained space distri¬ 
bution is not very sensitive to the value of a. 


normal kernel method (adopting = 0.30r, cf. Section 3.1.1); 
the spatial distribution of OB runaways producing bow-shocks 
(top panel Fig. is divided by the spatial distribution of the 
whole sample of OB runaways, including the HMXBs (bot¬ 
tom panel Fig. ||). Some OB runaways can be expected to not 
produce a bow-shock, because of their high space velocity or 
their large separation from the galactic plane. For that reason, 
the objects used for determining the probability distribution in 
Fig. are only those with intermediate velocities (5 km s ^ 
< Kpace < 80 km s“^), and those that are located within 
300 pc of the galactic plane. Note that areas with a small num¬ 
ber of runaways will have a rather inaccurate value of the bow- 
shock probability density (division through a small number). 
Fig. 1^ illustrates that the obtained distribution is not very sen¬ 
sitive to the value of a; for = 0.5r the space distribution 
becomes uniform at distances beyond ~ 0.7 kpc from the Sun. 

It can be seen that the distribution of nearby hot bubbles 
(from H 1 measurements) is consistent with the distribution of 
regions showing a low fraction of bow-shock producing run¬ 
aways. At larger distances this method of mapping hot regions 
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becomes inaccurate due to the rapidly decreasing number of 
objects. For example, the probability density maximum in the 
lower left corner of Fig. (top panel) is completely deter¬ 
mined by only one object beyond 3 kpc in that quadrant, an 
object that happens to show excess infrared emission. 

5.2. HMXBs have a lower probability to produce bow-shocks 
than dynamically ejected OB runaways 

Only one of the 11 FlMXBs has an associated wind bow-shock. 
Though small-number statistics, does this indicate that HMXBs 
have a lower probability to form a wind bow-shock? The bow- 
shock detection rate for Van Buren’s sample is about 40 %. If 
we take into account that the detection sensitivity of IRAS is 
poor for objects beyond 2 kpc, then we are left with 1 out of 5 
HMXBs. These numbers do not provide statistically significant 
results, but the different evolutionary origin, and therefore on 
average different kinematical age, of HMXBs and dynamically 
ejected runaways might result in a difference in detection rate 
of wind bow-shocks. 

Massive stars ejected by dynamical interactions should 
have escaped the OB association at a relatively early stage, 
when the cluster is still dense and the probability for close en¬ 
counters high. On the other hand, OB stars in HMXBs went 
through several stages of binary evolution before they became 
runaways after the occurence of a supernova within the sys¬ 
tem. Therefore, the kinematical age (i.e. the travel time since 
ejection took place) of HMXBs should be smaller than those 
of dynamically ejected runaways. As a consequence, HMXBs 
have a higher probability to be still enclosed in the hot and 
rarefied regions (e.g. superbubbles) that surround OB associ¬ 
ations. Further study is required to prove this hypothesis with 
observations. 

6. Summary and conclusions 

Based on our study of OB runaways and the phenomenon of 
wind-bow-shock formation we conclude that about 40 % of the 
OB-runaway stars produce wind bow-shocks. This fraction is 
a little bit higher than proposed by Van Buren and coworkers, 
but takes into account the estimated reduction in IRAS sensi¬ 
tivity required to detect a bow-shock. We investigate why the 
majority of OB runaways do not produce a bow-shock. The 
spectral type (temperature, luminosity, mass-loss rate), space 
velocity and separation from the galactic plane seem to be of 
minor importance. We propose that the temperature and density 
of the ambient medium are the dominant factors. In hot bubbles 
the speed of sound is larger than the typical velocity of a run¬ 
away star (~50 km s“^), so that runaways in these regions do 
not move supersonically and would not produce a bow-shock. 
OB-runaways produced through the binary supernova scenario 
should on average have a shorter kinematical age than those 
produced through the cluster ejection scenario, and thus have a 
higher probability to be still inside a hot bubble. For some run¬ 
aways, e.g. Cyg X-1, there is observational evidence that the 
system indeed is contained in a (super)bubble. 
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